Accepted for Publication in the Astrophysical Journal 

Preprint typeset using MTgX style emulateapj v. 6/22/04 



o 
o 



THE ASSEMBLY OF DIVERSITY IN THE MORPHOLOGIES AND STELLAR POPULATIONS 

OF HIGH-REDSHIFT GALAXIES ^ 

Casey Papovich 

Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721; papovich@as.arizona.edu 

Mark Dickinson^'^ 

National Optical Astronomy Observatory, 950 North Cherry Avenue, Tucson, AZ 85719; med@noao.edu 

Mauro Giavalisco 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; mauro@stsci.edu 

Christopher J. Conselice^ 

i/^ i Palomar Observatory, California Institute of Technology, MS 105-24, Pasadena, CA 91125; cc@astro.caltech.edu 

AND 

Henry C. Ferguson^ 

space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; ferguson@stsci.edu 
^ , Accepted for Publication in the Astrophysical Journal 

d 

ABSTRACT 

. We have used deep images from the Hubble Space Telescope to study the evolution in the morphologies, 

sizes, stellar-masses, colors, and internal color dispersion of galaxies in the Hubble Deep Field North. We 
selected two galaxy samples at 0.7 < z < 1.4 and L9 < z < 3.0 from a near-infrared, flux-limited catalog 
with complete photometric and spectroscopic redshift information. At z ~ 1 the majority of galaxies with 
M(B) < -20.5 have rest-frame optical morphologies of early-to-mid-type Hubble-sequence galaxies, and 
\ many galaxies show strong transformations between their rest-frame UV and optical morphologies. In stark 

contrast, galaxies at z ~ 2.3 all have compact and irregular rest-frame optical morphologies with little differ- 

■ ence between their rest-frame UV-optical morphologies, and with no clearly evident Hubble-sequence candi- 
i/^ \ dates. The mean galaxy size increases from z ^ 2.3 to 1 by roughly 40%, and the number density of galaxies 

■ larger than 3 kpc increases by a factor of ^ 7. The size-luminosity distribution of z ~ 1 galaxies is broadly 
consistent with that in the local universe, allowing for passive evolution. However, we argue that galaxies at 
z ~ 2.3 are not the fully formed progenitors of present-day galaxies, and they must continue to grow in both 
their sizes and stellar masses. We have quantified the differences in morphology by measuring the galaxies' 

O ' internal UV-optical color dispersion, which constrains the amount of current star-formation relative to older 

stellar populations. The mean and scatter in the galaxies' UV-optical total colors and internal color dispersion 
c/) ■ increase from z~2.3to 1. Atz~ 1 many galaxies with large internal color dispersion are spirals, with a few 
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irregular and interacting systems. Few z ~ 2.3 galaxies have high internal color dispersion, and we infer that 
those that do are also actively merging. Using simple models, we interpret the change in the total color and 
internal color dispersion as evidence for the presence of older and more diverse stellar populations at z ~ 1 that 
are not generally present at z ^ 2. The z ~ 2.3 galaxies do not increase their stellar diversity as rapidly as they 
could given basic timescale arguments and simple models. We conclude that the star-formation histories of 
galaxies at z ^ 2 are dominated by discrete, recurrent starbursts, which quickly homogenize the galaxies' stellar 
content and are possibly associated with mergers. The increase in the diversification of stellar populations by 
z ~ 1 implies that merger-induced starbursts occur less frequently than at higher redshifts, and more quies- 
cent modes of star-formation become the dominant mechanism. This transition in the mode of star formation 
coincides with the emergence of Hubble-sequence galaxies, which seems to occur around z ~ 1.4. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: stellar content — galaxies: structure 



1. INTRODUCTION 

Galaxy morphologies stem from complex formation and 
evolutionary records. The fact that present-day galaxies span 
such a great range of morphological type and that galaxy mor- 

^ Based on observations taken with the NASA/ESA Hubble Space Tele- 
scope, which is operated by the Association of Universities for Research in 
Astronomy, Inc. (AURA) under NASA contract NAS5-26555. These obser- 
vations are associated with programs GO/DD-6337 and GO-7817. 

^ also Space Telescope Science Institute 

^ also Department of Physics and Astronomy; The Johns Hopkins Univer- 
sity, 3400 N. Charles St., Baltimore, MD 21218 
4 NSF Fellow 



phology c orrelates with optica l luminosity, color, size, and 
mass re.g..lRoberts & HavneJl994.) indicates that galaxy ap- 
pearances depend strongly on their assembly histories. The 
quantitative study of the morphological properties of galax- 
ies at high redshifts has advanced rapidly due primarily to 
the high-angular resolution of the Hubble Space Telescope 
(HST). The HST instruments provide imaging with angular 
resolution of ~ Of'l, which corresponds to physical scales of 
less than 1 kpc for cosmologically distant galaxies. 

Many surveys using HST imaging have broadly 
shown that while the local Hubble Sequence se e ms to 
exist for z ^ 1-1.5 (e.g., Ivan den Bergh et all Il996t 
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Ravindranath et al. 2004), the morphological mix of 
galaxies at these redshifts shows an increase in the abundance 
of irregular and pec uliar galaxies, particularly at fainter mag- 
nitudes re.g.lPnver . Windhorst, & Griffiths 1995; ElHs 1997; 
iBrinchmann etalJ 11998: ,Brinchmann & EUis, ,200a) . How- 
ever, the formation and evolution of such systems is unclear. 
Several studies have concluded that a large fraction of giant 
spiral and other ea rly-t ype galaxies w ere alre ady in place by 
moder ate redshifts Jvan den Bergh et all 1 996t I Abraham et al. 
199a iBri n chman n'etal.i .1998:" ,Lillv et alJ 11998: Im et al. 
1999|: ISim ard et al J Il999t llm et al.l l2002t iM oustakas et al. 
2004; Ravindranath et al. 2004; Stanford et al. 2004) and that 
the characteristic galaxy size has changed little since z ^ I 
iLillyetal. 1998). Moreover, little discernible evolution 
is evident in such systems even when uncertainties in the 
morphological classification and oth e r selection effects 
are included ( Marleau & Simard 1998: Sima rd etalJ 11999: 
iRavindranath et al. 2004), which reinforces the conclusion 
that some fraction of the Hubble-sequence galaxy population 
had assembled by z ~ 1 . 

How these large, symmetric Hubble-sequence galaxies 
form and assemble their diverse stellar populations remains 
enigmatic. Logically, they descend from higher-redshift pro- 
genitors, but it is difficult to link the populations. The Hub- 
ble sequence seems to no longer apply at z ^ 1.5-2 (van 
den Bergh et al. 2001; van den Bergh 2002). The observed 
morphologies of galaxies at z ^ 2 are generally highly ir- 
regular or centrally compact (Giavalisco, Steidel, & Mac- 
chetto 1996; Lowenthal et al. 1997; Dickinson 2000; Daddi et 
al. 2004), bu t some large di sk-galaxy can didates exist (e.g., 
iLabbe etal .1 12003: Conselic e'et al.ll2004bh . The sizes, spa- 
tial densities, and stellar masses of galaxies at these redshifts 
are comparable to those of present-day bulges an d low er- 
redshift gal actic spheroidal components ( Giavalisco et al. 
, 1998t iLo wenthal et al.'ri9 97: .Giav alisco & Dickinson 
j iPapovich . Dickinson. & F ergusonl l2001: Shapley et al. 
2001^, and thus it is tempting to equate them with a 
one-to-one correspondence. However, the space den- 
sity of galaxies at z ~ 2 - 3 may be inconsistent with 
the fraction of bulge-dominated galaxies observed at mod- 
erate redshift (z ~ 1) in the HDF-N ( Marleau & Simard 
ri998h . Moreover, although the stellar-mass density at 
z ^ 1 is dominated by red, evolved galaxies (Bell etal. 



120041) . the majority of the galaxy stellar mass at z ^, 
2 exists in blue, irregular galaxies (e.g., Dickin son et al.l 
|2003; Conselice. Blackburne. & Papovich 2004a) with some 
contribution from dust-enshrouded, massive starbursts 
fporster-Schreiber et al. 2004; van Dokkum et al. 2004). 

The physical processes that govern how galaxies observed 
at z ^ 2 evolve into galaxies at more modest redshifts (z ~ 1) 
are poorly understood. This is partly due to difficulties in 
comparing the properties of galaxies observed at different 
rest-frame wavelengths and with different telescope and de- 
tector combinations. It is entirely expected that all observable 
galaxy properties (e.g., luminosities, colors, masses, mor- 
phologies, densities) evolve substantially over galaxies' life- 
times, and the subsequent observational biases further frus- 
trate studies of galaxy evolution when examining galaxy pop- 
ulations based on some selection criteria. 

One complication in the studies of galaxy morphology re- 
sults from the fact that a fixed bandpass samples lower rest- 



frame wavelengths with increasing redshift. For example, ob- 
servations at optical wavelengths probe only the rest-frame 
UV portions of the spectral energy distributions (SEDs) of 
galaxies at z ^ 1, and there is little relation between rest- 
frame UV and optical properties of galaxies. Local galaxies 
have been shown to exhibit strong transformations between 
their mor phologies at rest-frame UV and optical wavelengths 
(see, e.g. Kuchinski et al ^^200^, '2OOP; 'M arcum etal1l200TI: 
Windhorst et al. 2002; Papovich et al. 2003). The UV emis- 
sion from galaxies is generally dominated by the light from 
young, massive stars, and as a result a galaxy's morphology 
at UV wavelengths can appear as a later-type system com- 
pared to that derived at optical wavelengths. The composition 
(ratio of the number of young to old stars) and configuration 
(heterogeneity in the distribution of young and old stars) of 
the stellar populations within elliptical, spiral, and irregular 
galaxies varies widely. Galaxy morphology can vary strongly 
from optical to UV wavelengths, and comparing the morpho- 
logical properties of galaxies at UV wavelengths likely pro- 
vides an incomplete evolutionary picture. 

To compare the appearances between high-redshift and 
local galaxies, one requires observations at near-infrared 
(NIR) wavelengths, which probe rest-frame optical wave- 
lengths for z ^ 1 galaxies. Several recent studies have 
focused on the morphological properties of high-redshift 
galaxies in their rest-frame optical using HST high an- 
gular reso lution NICMOS d ata (e g. Teplitz et al. 1998; 
iDickinsonI 120001: "Corbin et al." *200U Moth & Elstoi; .2002: 
Consel iceeta]j 2003 : Stanford et al. 2004) and ground-based 
data (e.g. Labbe et al. 2003; Truiillo et al. 2004), and this was 
a strong motivation behind our NICMOS imaging survey of 
the HDF-N. This dataset combined with the HST/WFFC2 
data (Williams et al. 1996) provides high angular resolution 
images (~ O.^'l) in six bandpasses spanning 0.3-1.6 /im, 
which permits comparisons between galaxy morphologies at 
rest-frame UV-optical wavelengths as a function of redshift 
forO<z< 3. 

In this paper, we investigate the distribution of the mor- 
phologies, sizes, colors, and internal colors of high-redshift 
galaxies (z ~ 0.7 - 3.0) in the HDF-N using data available 
from the HST datasets. In a previous paper (Papovich et al.| 
2003, hereafter P03), we presented a diagnostic to mea- 
sure a galaxy's internal color dispersion, which quantifies 
the differences in galaxy morphology as a function of ob- 
served wavelength (i.e., it quantifies the morphological K- 
correction). Here, we apply this statistic to the rest-frame 
UV-optical colors of a sample of galaxies selected from the 
deep HST/mCMOS observations of the HDF-N. In § 2, 
we present a summary of the data and reduction procedures 
for the HDF-N data and define galaxy samples used here. 
In § 3, we describe the galaxy morphologies, and quantify 
their sizes, colors, and internal color dispersion. In § 4, 
we discuss differences in these quantities versus redshift, 
and compare these with simple scenarios of galaxy evolu- 
tion. In § 5, we summarize our results. Throughout this 
work, we assume a cosmology with Qm = 0.3, A = 0.7, and 
a Hubble parameter of /z = 0.7 (where the present-day Hub- 
ble constant is //q = 100/z k m s~^ Mpc~^). All magnitude s 
have units in the AB system (lOke. Gunn. & Hoess 'ell ll996h . 
mab = 3 1 .4 - 2.5 log((/jy )/ 1 nJy). We also refer to the four 
WFPC2 (F300W, F450W, F606W, F814W) and the two NIC- 
MOS (FllOW, F160W) bandpasses used on the HDF-N as 
^300, ^450, ^606, hi4, ^110, ^160, respectively. 
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2. THE DATA AND SAMPLE DEFINITION 

The complete HDF-N was observed with NICMOS for GO 
proposal 7817 during the UT 1998 June NICMOS refocus 
campaign.^ The details of the observations, data reduction, 
and object cataloging will be presented elsewhere (see also 
Dickinson et al. 2000; Papovich et al. 2001; Dickinson et al. 
2003; Stanford et al. 2004). Briefly, we observed the com- 
plete HST/WFFC2 field of the HDF-N with eight sub-fields 
with the NICMOS Camera 3. Each sub-field was imaged 
three times through both the FllOW and F160W filters, and 
each image was dithered into three separate exposures (result- 
ing in nine independent, dithered positions in each sub-field), 
with a net exposure time of 12600s per filter. The data were 
processed using standard STScI pipeline routines and custom 
software, and combined into a single mosaic that was regis- 
tered to the HDF-N WFPC2 images (WilHams et al. 1996) us- 
ing the "drizzling" method of Fructher & Hook (2002). The 
WFPC2 and /no images were then convolved to match the 
//i6o point-spread function (FWHM ^ 0'/22) in order to en- 
sure that photometry was measured through similar apertures 
in each image. Our tests on point sources in the PSF-matched 
images indicate that photometry between bands is accurate to 
within 5% for aperture radii greater than O.^^l. O bject cat- 
alogs were constructed using SExtractor iBertin & Arnouts 
Il996h by detecting objects in a combined /iio+/^i6o image 
and measuring fluxes through matched apertures in all in- 
dividual bands. As in Stanford et al. (2004), we will refer 
to the complete galaxy catalog as the HDF NICMOS Mo- 
saic, or HNM, catalog. Relative photometry from the HST 
WFPC2 and NICMOS data was measured through matched 
elliptical apertures defined using the radial moments of each 
galaxies' FllOW + F160W light profile (i.e., the SExtractor 
MAG_AUTO flux measurements). The difference between 
the flux measured in the total (MAG_AUTO) apertures and 
the isophotal (MAG_ISO) apertures has a median value of 
^ 6% for all of the galaxies with Hi^o <25.5. 

In this work, we define two galaxy samples from the HDF- 
N data selected on the basis of redshift and limiting magni- 
tude: a moderate-redshift sample, 0.7 <z< 1.4 (with a mean 
redshift z — 1); and a high-redshift sample, 1.9 < z < 3.0 
(with a mean redshift z — 2.3). We have used the most com- 
plete collection of HDF-N spectroscopic redshifts known (see 
Dickinson et al. 2003; and references therein). For the re- 
maining galaxies, we have used the photometric-redshift cat- 
alog of Buda variet al.> (2000) for the HDF-N based on the 
combined WFPC2 and NICMOS data. 

Because we analyze the two-dimensional galaxy appear- 
ances, we have selected galaxies on the basis of signal-to- 
noise ratio (S/N) down to a limiting isophotal level, which is 
more appropriate than a pure, flux-limited sample as galax- 
ies of larger sizes have lower S/N for a fixed flux density. 
We also select the objects in the bandpass that approximately 
corresponds to the rest-frame 5-band emission. We select 
all galaxies at z ~ 1 with S/N(/iio) > 20 within a limiting 
isophote above fi(Jno) ^ 26.5 mag arcsec"^. Similarly, we 

^ During the 1998 refocus campaign (HST Cycle 7N) the HST secondary 
was moved to the optimal focus for the NICMOS camera 3. This has not been 
done in cycles following the installation of the NICMOS Cooling System 
(NCS) during the HST servicing mission 3B. As a result, the image quality 
of the HDF-N mosaic is significantly improved over what is currently avail- 
able with NICMOS post-NCS. For example, the PSF FWHM of the HDF-N 
mosaic, ~ 0^/22 in Hi^q, is approximately 60% smaller than that of similar 
observations of the Hubble Ultra-Deep Field (R. L Thompson, private com- 
munication; see also Roye et al. 2003). 



select galaxies at z ~ 2.3 with S/N(//i6o) > 20 within a lim- 
iting isophote above /i(ifi6o) < 26.5 mag arcsec"^. Our tests 
have shown that these limits on object S/N within an isopho- 
tal aperture are reasonable for measuring galaxy morpholog- 
ical and structural parameters (see P03 and discussion be- 
low). These limits in S/N correspond to approximate mag- 
nitude limits measured within the elliptical apertures (i.e., the 
SExtractor MAG_AUTO) of /no - Hm - 25.5. The final 
galaxy samples have 113 objects at z ~ 1 and 53 objects at 
z ^ 23 that satisfy these selection criteria in S/N and red- 
shift. Of the samples, 69% of the z ~ 1 galaxies and 62% of 
the z ~ 2.3 galaxies have spectroscopic redshifts. In the fol- 
lowing analysis, we also compare the properties of the galaxy 
samples down to a fixed absolute magnitude, M(B) < -20.0, 
where both samples are mostly complete (see Dickinson et al. 
2003). This luminous sub-sample contains 41 objects at z ~ 1 
and 46 objects at z ~ 2.3. 

The redshift ranges used here were chosen to take advan- 
tage of two attributes of the data. The apparent galaxy size re- 
mains roughly constant for 0.7 < z ^ 3.0 because the angular- 
diameter distance — the ratio of an object's physical size to 
its angular size — changes little (< 18%) for the default cos- 
mology: from 7.1 kpc arcsec"^ at z = 0.7 to 8.5 kpc arcsec"^ 
at z ~ 1 .6 (this is generally true for any plausible cosmology). 
Thus, the HST images provide roughly constant angular res- 
olution of ~ 1.5 kpc for all the galaxies in this dataset nearly 
independent of redshift. Note, however, that the co-moving 
volume contained in the redshift slice 1.9 < z < 3.0 is 2.3 
times than contained by 0.7 < z < 1.4. 

The other advantage comes from the fact that at these red- 
shifts combinations of the WFPC2 and NICMOS bandpasses 
probe approximately the same wavelengths in the rest-frame, 
which allows us to minimize morphological ^-corrections 
between the two samples. We broadly define two rest-frame 
wavelength regions: the mid-UV (~ 2400-3800 A), and op- 
tical (- 4000 - 6000 A). For the 0.7 < z < 1 .4 galaxy sample, 
the mid-UV-, and optical-wavelength regions correspond to 
the V606 and /no bandpasses, whereas at 1.9 < z < 3.0 these 
rest-frame wavelength regions correspond to the and Hieo 
bandpasses. We note that using a single bandpass to track 
a specific rest-frame wavelength for a wide redshift interval 
does induce overlap in wavelength and redshift coverage due 
to the width and central wavelengths of the of the HST filter 
curves. For example, due to its large width, the centroid of 
the /no filter samples Ao - 4600-6500 A for z 0.7-1.4, 
and Ao - 2750-3800 A at z - 1.9-3.0. One could devise 
a method to interpolate the observed bandpasses to acquire 
the flux at a common wavelength. However, this would in 
principle hinder the analysis here, as one of the goals of this 
work is to contrast and compare the individual color features 
of galaxies for common rest-frame wavelengths. Interpola- 
tion between bands would only suppress any intrinsic signal. 
Therefore, we have opted to use the fixed bandpasses as ap- 
proximate surrogates for the rest-frame colors. Furthermore, 
our tests on the UV-optical colors of local galaxies (P03) has 
shown that the internal color dispersion is most sensitive to 
changes due to the the heterogeneity of the stars that con- 
tribute to the flux redward and blueward of the B aimer and 
4000 A breaks. For the galaxies considered here, these breaks 
are always confined between the chosen bandpasses. 
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3. GALAXY MORPHOLOGICAL AND STRUCTURAL PROPERTIES 

Figures [I] and |2| show images of the 28 most luminous 
galaxies (in the rest-frame 5-band) in each of the z ~ 1 and 
z ~ 2.3 samples. These figures illustrate the morphologies 
of the brightest galaxy members of each sample. Although 
the galaxies in each sample span comparable rest-frame lu- 
minosities, -20.5 > M(B) > -22, their morphologies are dis- 
tinctly different. Many of the galaxies at z ~ 1 have regular 
and symmetric morphologies, particularly in the rest- frame 
optical images, and are similar to "normal" Hubble-sequence 
galaxies. Our by-eye classification of the /no morpholo- 
gies finds that of the brightest z ~ 1 galaxies 0:^43% (12/28) 
have elliptical, lenticular, or spheroidal morphologies (HNM 
144, 789, 847, 886, 1014, 1091, 1031, 1214, 1356, 1418, 
1453, 1523), 43% (12/28) have spiral or disk morpholo- 
gies (HNM 214, 388, 448, 577, 779, 1120, 1168, 1213, 1240, 
1348, 1488, 1525), and only 14% (4/28) have irregular or 
compact morphologies (HNM 333, 360, 909, 1316). Many of 
the objects show strong differences between their morpholo- 
gies at UV and optical wavelengths (i.e., strong morpholog- 
ical ^-corrections), which is particularly apparent in those 
galaxies with spiral and/or disk morphologies. 

In contrast, the brightest galaxies at z ~ 2.3 have either ir- 
regular (HNM 552, 1357, 163, 1358, 813, 814, 1513, 109, 
850, 758, 701, 502, 804) or compact morphologies (HNM 1 10 
661, 843, 466, 274, 741, 62, 67, 1047, 1550, 503, 229, 272, 
230, 1 199). Unlike the galaxies at z ~ 1, there are no galaxies 
with disk-like or spiral-galaxy morphologies at z ~ 2.3: no 
galaxies show evidence for spiral arms or for being edge-on 
disk systems. This is broadly similar to the morphological 
descriptions o f z ^ 3 galaxies based solely on rest-frame UV 
observations JGiav alisco et al. 1996; Lowenthal et al. 1997). 
The majority of the optically bright galaxies at z ~ 2.3 have 
similar morphologies in the V606- (rest-frame 1500 A) and 
//i6o- (rest-frame 4400 A), which is even more significant 
because this wavelength baseline is longer than that shown 
for the z ~ 1 sample in Figure [I] The lack of significant 
morphological ^-correction combined with the fact that these 
galaxies have blue UV-optical colors (Dickinson et al. 2003) 
implies that recent star formation dominates the UV-optical 
morphologies, which we discuss in more detail below (see 
also the discussion in P03). 

3.1. The Galaxy Size-Luminosity Relationship 

The panels in figures [l] and |2| show the images of the z ~ 1 
and 2.3 galaxies at the same physical scale (i.e., each panel 
is approximately 40 kpc x 40 kpc). Qualitatively, the typical 
galaxy size appears to increase from z~2.3toz~l, with the 
most dramatic evolution in the rest-frame 5-band images. In 
this section, we quantify this observed evolution in terms of 
the galaxy size-luminosity distribution. 

We derived approximate absolute 5-band magnitudes by 
interpolating the V6064 14-^110^160 photometry to measure the 
apparent magnitude at rest-frame 4400 A and applying the 
distance moduli with the default cosmology for the observed 
redshift. We quantify galaxy sizes with the half-light radius 
ri/2, which is defined as the radius of the aperture that encom- 
passes 50% of the object's total flux as derived by the SExtrac- 
tor software. Our tests using simulated data with noise proper- 
ties similar to that of the HDF-N indicates that the SExtractor 
measurements of the half-light radius are accurate (< 10% 
rms) for S/N > 20, which is satisfied for the galaxy samples 
used here. 



Figure [3l shows the distribution of the derived half-light 
radii for the HDF-N galaxies as a function of their ab- 
solute magnitudes at rest-frame 4400 A. The upper en- 
velope of rest-frame B-band luminosities is roughly the 
same at z ^ 1 and z ^ 2.3 (see also the discussion above; 
iDickinso n et al jl2003h . These luminosities are comparable to 
those of the larg est present-day early-type spirals and ellipti- 
cals (^e Vaucou leurs et alJri99U iBender. B urnstein. & FabeJ 
1992). Many galaxies in both samples have luminosities well 
in excess of the present-day characteristic "L*",^ although the 
majority of galaxies extend to fainter luminosities. 

The average galaxy half-light radius grows with time from 
z ~ 2.3 to 1. The z ~ 2.3 galaxies have relatively small in- 
trinsic sizes, 1 ^ ^1/2 ^ 4 kpc, with a mean value {ri/2) — 
2.3 ± 0.3 kpc for M(B) < -20.0, where the error represents 
the uncertainty on the mean. In contrast, the measured radii 
of the z ~ 1 galaxies are larger for comparable absolute mag- 
nitudes. The mean half-light radius for galaxies with M(B) < 
-20.0 is (ri/2) — 3.2 ± 0.3 kpc, which is roughly 40% higher 
than the average size of the z ~ 2.3 galaxies. We find that 
based on the two-dimensional Kol mogoro-Smirnov (KS) test 
JPeacockll983l:lFasano & Francesc hini 1987) the z ~ 2.3 and 
z ~ 1 galaxy samples, for M(B) < M*(B), trace different dis- 
tributions at the 99.96% confidence level (c^ 3.6a). The 
luminosity-size distribution is qualitatively similar to results 
reported by Trujillo et al. ( 2004) for the HDF-S, but a quanti- 
tative comparison is beyond the scope of this work given the 
differences in rest-frame bandpass and sample redshifts. 

In addition, there is a dramatic increase in the number den- 
sity of large galaxies from z ~ 2.3 to 1. Considering only 
galaxies with ri/2 > 3 kpc and M(B) < -20.0, we find 24 
galaxies at z ~ 1 and 8 at z ~ 2.3. Because the co-moving 
volume of the z ~ 2.3 redshift range is larger by a factor of 
2.3, there is an increase in the number of large galaxies by a 
factor of ~ 7 between these redshifts. Increasing the size limit 
to > 4 kpc makes the evolution even more dramatic: there are 
13 such galaxies at z ~ 1 compared with only 2 at z ~ 2.3, but 
this comparison suffers from obvious small-number statistics. 
Note that our magnitude limits are also conservative. If we in- 
clude fainter sources (e.g., to M(B) < -19), then the observed 
change is even greater: the number density of sources with 
ri/2 > 3 kpc increases by a factor of 11 from z ~ 2.3 to 
z ~ 1, but at the fainter magnitudes the z ~ 2.3 sample proba- 
bly suffers from some incompleteness. 

Note that we made no attempt to remove the contribution 
of the image PSF from the measured half-light radii. The 
NIC3 F160W PSF has a small intrinsic half-light radius, 
^i/2,PSF — 0. 15^', which for these redshifts and default cosmol- 
ogy corresponds to ri/2 — 1.2-1.4 kpc. The majority of ob- 
jects have ri/2 ^ 2 kpc, to which the PSF contributes < 20%. 
However, the PSF dominates the extrinsic sizes for objects 
with ri/2 ^ ^i/2,PSF- Because the PSF can only increase the 
measured half-light radii, its effects will tend to decrease the 
amount of intrinsic size evolution between z ~ 2.3 and z ~ 1. 
Therefore, if anything our result underestimates the intrinsic 
evolution. 

There is a deficit of HDF-N galaxies at 0.7 < z < 1.4 in 
the HDF-N with ri/2 ^ 8 kpc. Although this result could 
be interesting for the evolution and growth of galaxies for 
z ~ - 1 , it likely stems from the fact that the HDF-N data 
probes a small volume. The present-day bivariate (radius- 

^ Norberg et al. 2002 find a present-day characteristic luminosity of 
M{BY = -20.5 [AB], after slight adjustments in cosmology and bandpass. 
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Fig. L — Rest-frame UV and optical images of 28 galaxies in the 0.7 < z < 1.4 sample with the highest luminosities in the rest-frame 5-band from the 
HDF-N. The left panel for each galaxy shows the WFPC2 F606W image, which corresponds to the rest-frame mid-UV wavelengths. The right panel shows 
the NICMOS FllOW image, which depicts each galaxy approximately in the rest-frame 5-band. Galaxies span luminosities -22.7 < M{E) < -20.4, and have 
redshifts as indicated in each image. Redshifts shown within parentheses indicate cases where only photometrically derived redshifts are available. Note the 
diversity in the galaxies appearances as a function of rest-frame wavelength. Each panel is 5 x 5 arcsec^ , which corresponds to ~ 40 x 40 kpc^ in the rest-frame. 



, luminosity-) distribution function fde Jong & Lacevll20QCl) 
predicts that there should be 1.6 ± 1.3 objects with ri/2 ^ 
8 kpc in the co-moving volume of the HDF-N for 0.7 < z < 
1.4. Therefore, the lack of galaxies with large radial sizes 
in the HDF-N at these redshifts is not strongly inconsistent 
with the present-day galaxy population: the volume is sim- 
ply too small to expect many giant galaxies. This conclusion 
is reinforced by th e presence of larg e galaxies in th e disk- 
galaxy sample s of iLilly et al J (Il998h . Isimard etal.l (|1999'), 
andlRavindran ath et al.l (120041) to z ^ 0.7- 1, which are drawn 
from redshift surveys that cover larger solid angles. However, 
if the local galaxy size-luminosity distribution were applica- 
ble at higher redshifts, then one would expect to find 5 ± 2 
galaxies with 1 .9 < z < 3.0 and ri > 8 kpc, and none are ob- 
served. Note that although Labbe et al. (2003) report several 
galaxies with ri/2 — 4-6 kpc at z ~ 2-3 in the HDF-S, there 
are no candidates at these redshifts with larger radii. 

The (rest-frame) 5-band half-light radii of the z ~ 2.3 
galaxies are consistent with those reported from z ~ 3 



galaxies from observations solely in the rest-frame UV 
jGiavalisco et al. 1996; Lowenthal et al. 1997). This implies 
that UV-selected galaxies at these redshifts have high surface 
brightness at both rest-frame UV and optical wavelengths, 
and are not generally UV-bright star-forming regions within 
otherwise normal, larger host galaxies. 

3.2. Effects of Surface Brightness on Galaxy Sizes 

We have considered the effects of surface-brightness dim- 
ming, which rapidly increase with redshift (the bolometric 
surface brightness scales as (1 +z)~^). We simulated images of 
the luminous (M(B) < -20.0) galaxies from the z ~ 1 sample 
as they would appear at z = 2.7 (i.e., somewhat higher than 
the mean redshift z = 2.3 in order to achieve more conserva- 
tive limits). The general details of this process are presented 
in P03, but briefly, we resampled the NICMOS /no images 
of the z ~ 1 galaxies to the physical pixel scale for z = 2.7. 
We then appropriately reduced the pixel intensity to account 
for surface-brightness dimming and inserted them into blank 
regions in the NICMOS Hi^o images. Because the relative 
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Fig. 2. — Same as Figure[T] but here the figure shows the 28 galaxies from the 1.9 < z < 3.3 galaxy sample with the highest luminosities in rest-frame 
5-band. The left panel for each galaxy shows the WFPC2 F606W image, which corresponds to rest-frame far-UV wavelengths. The right panel shows the 
NICMOS F160W image, which corresponds approximately to the rest-frame 5-band (i.e., the same rest-frame band as FllOW for galaxies at z ~ 1). Note the 
lack of morphological transformation between wavelengths, even though these images span a longer wavelength baseline compared to Figure [l] The galaxies 
span luminosities -22.5 < M(B) < -20.7 (cf. to the similar range of galaxy 5-band luminosities shown in FigureQ see also Dickinson et al. 2003). The size of 
each panel is 5 x 5 arcsec^ , which corresponds to ~ 40 x 40 kpc^ in the rest-frame. 



resolution does not change substantially from z = 0.7 to 2.7, 
we have made the approximation that the image PSF does not 
vary in the simulated images. To compare these simulated 
z = 2.7 galaxies with the z ~ 2.3 HDF-N galaxy sample, we 
repeated the source detection and photometry processes using 
the SExtractor software to each simulated "redshifted" galaxy. 
The galaxy absolute magnitudes and characteristic radii were 
calculated on the simulated images in exactly the manner de- 
scribed previously. 

Figure|4]shows the distribution of galaxy sizes with M(B) < 
-20.0 for the galaxy samples, and for the the z ~ 1 sample 
simulated at z = 2.7. The majority of the luminous z ^ I 
galaxies would be detected if they were present at z ~ 2.7. 
The mean half-light radius of the z ~ 1 galaxies simulated 
at z = 2.7 is reduced slightly: (ri/2,sim) — 3.1 kpc (reduced 
by < 20% compared to the results above). Using the one- 
dimensional KS-test to the half-light distributions of the 
z ~ 2.3 sample and the z ~ 1 sample simulated at z = 2.7, 
we find that these populations stem from different parent dis- 



tributions at the 3.9(7 level. Similarly, applying the two- 
dimensional KS-test to the subset of bright galaxies with 
M(B) < -20 (as done previously), the radius-luminosity dis- 
tributions for the z ~ 2.3 galaxies and the z ~ 1 galaxies sim- 
ulated at z = 2.7 trace different underlying distributions at the 
99.6% confidence level (c^ 2.9a). Thus the strong size evo- 
lution that exists between galaxies at z ~ 2.3 and 1 is robust 
against surface-brightness dimming effects. This fact that the 
1 .9 < z < 3.0 redshift range contains more co-moving volume 
than that at 0.7 < z < 1.4 actually strengthens the significance 
of the result. For example, if we scale the number density of 
objects in each redshift interval to the same co-moving vol- 
ume, then the one-dimensional KS-test derives a > Sa like- 
lihood that the half-light distributions do not stem from the 
same parent sample. 

3.3. The Internal Color Dispersion of Galaxies 

In § 13.11 we describe the increase in the sizes of galaxies 
from z ~ 2.3 to 1 . If galaxies grow monotonically in size, then 
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Half-light Radius [kpc] 

Fig. 3. — Distribution of galaxy half-light radii versus absolute magnitude in rest-frame light at 4400 A. The absolute magnitudes are presented in AB 
units for the default cosmology. The plot shows the radius-luminosity distribution for the HDF-N z ~ 2.3 (open stars) and z ~ 1 (filled squares) galaxy 
samples, along with reference data collected from the RC3 catalog (de Vaucouleurs et al. 1991). The solid line indicates the relation for local E/SO galaxies 
teinggeli, Sandage, & Tarenghi 1984), and the dashed and dotted lines indicate the relation for local Sab/Sbc and Scd/Sdm/Irr galaxies, respectively IFreemanI 
[l97niiRoche et al.^1996 ). Note that the NICMOS F160W PSF corresponds to a half-light radius of 1.2- 1.4 kpc for z ~ 0.7-3.0. 



relatively young stars will form at larger radii at later times. 
An observable consequence of this growth may be that inter- 
nal color structure develops within galaxies at later epochs. 

To investigate possible internal color evolution in the HDF- 
N galaxy samples, we have quantified the galaxies' inter- 
nal color dispersion between (rest-frame) UV-optical wave- 
lengths of the HDF-N galaxies using the statistic defined in 
P03. The internal color dispersion quantifies the differences in 



galaxy morphology as a function of wavelength. In particular, 
the (rest-frame) UV-optical internal color dispersion charac- 
terizes diversification in the spatial configuration of young 
stellar populations (which dominate the rest- frame UV emis- 
sion) and older stars, which contribute more to the optical 
light. In summary (see P03 for details), the internal color dis- 
persion, ^, is defined as the dispersion of the pixel intensities 
between two bandpasses about a mean color, normalized to 
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Fig. 4. — Distribution of galaxy sizes for the z ~ 2.3 and z ~ 1 HDF-N samples for objects with M(B) < -20.0. In both panels, the shaded histograms show the 
measured distributions. In the right panel, the thick, solid lines illustrate the observed distribution of the z ~ 1 galaxies as measured from their appearances when 
simulated at z = 2.7 (see text), and normalized to the total number of objects in the shaded histogram of each panel. Note that the redshift range, 1.9 < z < 3.0, 
contains approximately 2.3 x the co-moving volume as 0.7 < z < 1.4, and so the deficit of large galaxies at z ~ 2.3 has even greater significance than what is 
illustrated in the figure. In each panel, the dotted line indicates the half-light radius of the NICMOS PSF. 



the total flux of the galaxy, i.e., 



E(h-f3)^-E(B2-aB,)^ 



(1) 



where h and h are the pixel intensity values from each image 
obtained in two bandpasses. The scaling factor a is the flux 
ratio, or color, between the images, while the linear offset 
P adjusts (if necessary) for differences in the relative back- 
ground levels of the two images. To account for internal color 
dispersion resulting purely from fluctuations in the sky back- 
ground, we subtract the contribution from background regions 
(i.e., image sections with no detected sources) in each image 
^1 and B2, respectively. The internal color dispersion, ^, is 
flux-independent, quantifies the morphological ^-correction 
of a galaxy between the bandpasses h and I2, and has an ana- 
lytically understood error distribution. 

We have computed the internal UV-optical color dispersion 
for the galaxies in the HDF-N samples using the V^oe and /no 
images for the z ~ 1 sample and hu and ifi6o images for the 
z ^ 23 sample. The HST images do not have a straightfor- 
ward distribution of background pixel values due to the cor- 
related pixels from the data-reduction process of "drizzling" 
the images, and due to the fact the the WFPC2 and NICMOS 
Fl low images were convolved to match the PSF of the NIC- 
MOS F160W image. Both properties of the data will tend 
to suppress inter-pixel noise. Thus, it is necessary to test that 
the analytic uncertainties in the internal color dispersion given 
in P03 accurately describe the statistical errors for galaxies 
in the HST data. We repeatedly inserted artificial galaxies 
with known internal color dispersion into blank regions of the 
HDF-N images, and then computed the distribution of mea- 
sured internal color dispersion values. The width of this dis- 
tribution provides an estimate of the uncertainty on ^. Indeed, 



our simulations indicated that the analytically derived error 
(5(0 had underestimated the intrinsic uncertainty. We subse- 
quently rescaled the derived 6(0 values by factors that were 
empirically derived from the simulations, i.e., S'(0 = C^(0, 
with C = 5.98 for aVeoeMo) and C = 7.14 for aku.Hml It 
should be noted that while the corrected errors provide a good 
estimate for the statistical errors arising from the background 
pixel-to-pixel noise, they likely still underestimate the real 
uncertainties as they neglect systematic effects arising from, 
e.g., inaccurate PSF matching, image registration, image cal- 
ibration, etc. In particular, errors in PSF-matching will have 
a stronger effect on compact, and centrally peaked objects, 
and are probably the dominant error source for the brightest 
objects. 

Figure|5lshows the distribution of internal UV-optical color 
dispersion as a function of redshift for the galaxies in the 
HDF-N samples. Large uncertainties on ^ are generally in- 
dicative of objects with low signal to noise in a particular 
bandpass. As a fiducial reference, we use the internal color 
dispersion values for two local galaxies observed in the same 
rest-frame bands, M81 (the early-type spiral galaxy with the 
highest internal color dispersion in the local P03 sample) and 
M74 (a "typical" mid-type spiral galaxy in the P03 sample). 
In Figure |3 we show the derived internal color dispersion 
for these two galaxies after convolving them with the NIC- 
MOS F160W PSF, and resampling them to the same physical 
pixel scale as the high-redshift HDF-N galaxies in the HST 
images. This illustrates the internal color dispersion in these 
fiducial local galaxies (with no evolution) relative to the high- 
redshift objects. 

There is an evident increase in the upper envelope of in- 
ternal color dispersion between z ^ 23 and 1. Formally, the 
mean and scatter increase by roughly a factor two from z ~ 2.3 
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Fig. 5. — The internal color dispersion measured between the (rest-frame) 
UV-optical colors for the HDF-N galaxy samples as a function of redshift. 
The plot shows the internal dispersion: ^(V^o6,J\\o) and C(4i4 5^i6o) for the 
z ~ 1 and z ~ 2.3 galaxies, respectively. Error bars indicate the approximate 
Icr uncertainties (see text). The dashed lines in each panel show the derived 
internal color dispersion for two local galaxies (M81 and M74) after convolv- 
ing with the NICMOS F160W PSF, and resampling to the same pixel scale 
as the high-redshift HDF-N galaxies. 



to 1: from {ahu.Hm)) ^ 0.02 at z - 2.3 to (e(y606,/iio)) ^ 
0.04 at z 1; and from cr(^[/8i4,i^i6o]) - 0.03 at z 2.3 
to cr(^ [/g 14, //i6o]) — 0.05 at z ~ 1. However, the change 
in the number of galaxies with large internal color disper- 
sion is more dramatic. Most of the galaxies at z ~ 2.3 have 
^(/8i4,^i6o) < 0.05 (46/53), with a small tail to ^ - 0.1 (7/53). 
In contrast, the z ~ 1 sample contains many more objects 
with high internal-color dispersion — 28/1 13 (25%) of galax- 
ies have ^(V606,'/iio) > 0.05 — including a sizable number 
(8/113) with values comparable to, or higher than, the mid- 
type spiral M74. 

Note that it is unlikely that we are missing some signifi- 
cant fraction of galaxies at z ~ 2.3 with high internal color 
dispersion. This is based on the facts that the 1.9 < z < 3.0 
redshift interval contains more than twice the co-moving vol- 
ume as that of the 0.7 < z < 1.4 (although cosmic variance 
uncertainties apply), the galaxies have been selected with the 
same limiting surface-brightness criteria, and the object reso- 
lution changes little for these redshifts (see § 2). Furthermore, 
in § 13. 41 we argue that surface-brightness dimming effects are 
not responsible for the change in the internal color dispersion 
at these redshifts, which reinforces our results from P03 that 
the intrinsic rest-frame UV-optical internal color dispersion 
of local galaxies is recoverable for < z ^ 3. This evolution 
is not dominated by surface-brightness-dimming effects. 

Figure shows the galaxies with the highest ^(Veoe^-^iio) 
in the z ~ 1 sample. Each object shows a significant quali- 
tative difference between the rest-frame UV and optical im- 
ages, which is quantified by the higher ^ values. Many of 
these galaxies have features indicative of early-to-late-type 
spiral galaxies (e.g., HNM 448, 860, 913, 1488, 1495), with 
a strong, central core in the /no images, and prominent star- 
forming spiral arms in the V606 images with relatively bluer 
colors. These objects have internal color dispersion and color 



residual images qualitatively similar to that observed in local 
spirals, e.g., M81, M74, MlOO (P03). Several of the z 1 
galaxies in the figure show the appearance of inclined disk 
systems (e.g., HNM 577, 685, 1076, 1447), and the resid- 
ual colors may suggest patchy dust and/or varied stellar con- 
tent similar to nearby analogs; e.g., M82 and UGC 06697 
(P03). There is also an actively star-forming galaxy with a 
disturbed morphology at all wavelengths and strong residual 
colors (HNM 909). This galaxy has a red core and blue outer 
regions, and it is possibly an spiral with some interacting ac- 
tivity, and a diverse stellar content and variable dust opacity. 
HNM 1453 shows clear multiple components of very different 
rest-frame UV-to-optical colors: a red elliptical galaxy with 
an off-center, blue component. This object has been identi- 
fied as a gravitational lens candidate of a background object 
(see, e.g., Zepf. Moustakas. & Davis 1997). Thus, it seems 
that this object consists of two objects along the line of sight, 
which were not split in the HNM catalog. We will exclude 
this object from further analysis. 

Figure El shows the galaxies in the z ~ 2.3 sample with the 
highest ^(/8i4,//i6o) values. Nearly all the galaxies show evi- 
dence for disturbed morphologies, and the internal color dis- 
persion probably arises from variations in the relative star- 
formation and variable dust obscuration. The exception is 
possibly HNM 163, which shows evidence for a redder core 
and a bright, off-center blue knot, and may result from varia- 
tions in the stellar content. However, this galaxy has evidence 
for high dust obscuration (Papovich et al. 2001), which may 
also dominate the internal color dispersion. 

Several pairs of the galaxies in z ~ 2.3 sample have been 
"deblended" by the SExtractor software into two components: 
HNM 162 + 163; HNM 813 + 814; HNM 1357 + 1358 
(see Figure |2||. These galaxies all have spectroscopic red- 
shifts, but the multiple galaxy components are merged at typ- 
ical ground-based resolution (^ I"), and as a consequence 
their spectra are blended as well. How ever, the p hotometric 
redshifts of the individual components ( Bu davari et alJI2000l) 
are approximately identical in each case (&/(l +z) ^ 0.02). 
Althou gh these individual components have similar redshifts 
(see Dickinson et al. 2003), it is unclear whether these objects 
actually are multiple components of the same parent galaxy 
or close galaxy pairs. In the analysis above we have kept 
the galaxies separated, but we have tested how merging them 
into single galaxies would affect the derived sizes and internal 
color dispersion. In three of the cases, the derived half-light 
radius increases by a factor less than two, with the exception 
being the case of HNM 162 + 163 where the half-light radius 
changes insignificantly. This would shift the galaxies in the 
radius-luminosity diagram (Figure|3ll, but would not alter any 
of the conclusions from that plot. In all four cases, the inter- 
nal color dispersion of the merged pairs equals essentially the 
luminosity-weighted average of the internal color dispersion 
of the individual components. This has little effect on the dis- 
tribution in Figure|5l and reduces the internal color dispersion 
of HNM 163 and 814, which have two of the largest values in 
the z ~ 2.3 sample. Using the values for the merged objects 
would actually augment the change in the scatter of the inter- 
nal color dispersion from z ~ 2.3 to 1. We therefore keep the 
objects separated in order to be conservative in this evolution. 

3.4. Effects of Surface Brightness on the 
Internal Color Dispersion 

In this section we consider the effects from surface- 
brightness dimming on the internal color dispersion for in- 
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Fig. 6. — The internal color differences of the HDF-N galaxies with the highest ^(V6065-^iio) values in the z ~ 1 sample. For each galaxy, the panels display the 
V6O6 and /no images (as labeled); the right-most panel shows the difference between the images after scaling to remove the mean total color, i.e., these residuals 
show I\ - ah -(3, where h , and h refer to the two images given in the panel inset as ^(/i 5/2), and a and (3 are the fitted constants (see text for full description). 
The corresponding internal color dispersion value is shown in each color-residual image. Note that the color residuals are not the same as the definition of 
internal color dispersion (the latter is the squared sum of the former; divided by the squared sum of the image flux). Light/Dark regions correspond to features 
with internal colors that are redder/bluer than the mean total color. Each panel has dimensions of 5 x 5 arcsec^ , which corresponds to 2:^ 40 x 40 kpc^ in the 
rest-frame for these galaxies. 



creasing redshift. In P03, we tested this for local galaxies 
observed by UIT, and concluded that the internal color dis- 
persion was fairly robust to z ^ 3 for galaxies with sufficient 
S/N in images with HST-XikQ angular resolution. Here, we 
expand these tests by simulating the appearances of bright 
(M(5) < -20.0) galaxies in the z 1 HDF-N sample. We 
simulated the appearance of each z ~ 1 galaxy at z = 2.3 (the 
median redshift of the high-redshift HDF-N sample) by ap- 
propriately resampling and dimming their images in the V6O6 
and /no bands, and adding noise from blank regions of the 
/814 and /fi6o images (see § 13.21 also P03). We then compute 
the internal color dispersion of these simulated galaxies using 
the procedures above. 

FigurelHlcompares the measured internal color dispersion of 
the z ~ 1 galaxies with the measured values from these galax- 



ies simulated at z = 2.3 in /gn and //160 images. In general, the 
internal color dispersion values of the simulated high-redshift 
galaxies correlate with their measured values at z ~ 1 . There 
is a slight reduction in the internal color dispersion values, 
which results from decreased signal-to-noise from cosmolog- 
ical surface brightness dimming. 

Quantitatively, the number of galaxies with internal color 
dispersion values above the fiducial value ^ > 0.05 drops from 
roughly 35% at z ~ 1, to 26% for the simulated galaxies at 
z = 2.3. However, this is still significantly greater than the 
15% of galaxies in the z ~ 2.3 sample with M{B) < -20 and 
^ > 0.05. Furthermore, the reduction of the internal color dis- 
persion in the simulated galaxies is possibly overestimated. 
In $ 14.21 we argue that the /814-//16O internal color dispersion 
at z ~ 2.3 is more sensitive to heterogeneities in a galaxy's 
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Fig. 7. — The internal color differences of the HDF-N galaxies with the highest C(/8i4,^i6o) values in the z ~ 2.3 sample. For each galaxy, the panels display 
the /gu and //160 images (as labeled); the right-most panel shows the difference between the images after scaling to remove the mean total color (see caption of 
fig|6|. Each panel has dimensions of 5 x 5 arcsec^, which corresponds to ~ 40 x 40 kpc^ in the rest-frame for these galaxies. 



Stellar populations compared to the V606-'^iio internal color 
dispersion at z ~ 1 . This is a consequence of the fact that at 
z = 2.3, the /814 bandpass probes rest-frame ^ 2400 A, while 
at z ~ 1 , the V6O6 bandpass is sensitive to slightly redder rest- 
frame wavelengths, ~ 3000 A. This is not taken into account 
in these simulations, and would increase the ^(/8i4,^i6o) val- 
ues of the simulated z = 2.3 galaxies by as much as ~ 40% 
Tsee^lOl. 

4. DISCUSSION 
4.1. The Growth of Galaxy Sizes 

The average size of luminous galaxies in the HDF-N be- 
tween z ~ 2.3 and z^l increases by roughly a factor of :^ 1 .4. 
Characterizing this as a simple power-law with redshift, we 
find that the galaxy sizes evolve as ri/2 ^ (1 for 
z ^ 1-2.3. Hierarchical models predict that the character- 
istic sizes of galaxy halos grow as, rvir(z) cx m\[^ H{zT^^^ cx 



VvivHizT^ re.g./Bouwens & Silk"2002VFerguson et a ll2004h . 
For the default cosmology, the Hubble constant evolves as 
H{z) ~ (1+z)^/^, which implies that galaxy sizes should 
evolve roughly as ry/r ~ (1 Assuming the half- 

light radius scales as the virial radius, the observations here 
support this hierarchical picture of matter assembly within 
halos for z ~ 1 - 2.3. This result broadly agrees with pre- 
dictions from other s emi-analytical hierarchical models of 
galaxy formation fe.g..'Baugh et al."1998':' Mo. Mao. & Whitd 
Il999; Somerville, Primack, & Faber 2001). 

The radius-luminosity distribution of z ~ 1 galaxies is 
roughly consistent with that of present-day galaxies and pas- 
sive luminosity evolution. For example, for a monotonically 
evolving stellar population model (e.g., using the models of 
Bruzual & Chariot 2003) formed in a burst of star-formation 
at z/ ~ 4, one expects AMg ^ 1 mag from z = 1 to z = 0. Thus, 
the data points for the z ~ 1 sample under this luminosity evo- 
lution would shift to regions of the Mb - ri /2 distribution that 
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Fig. 8. — The effects of surface-brightness dimming on the internal color 
dispersion. The figure compares the internal color dispersion values of galax- 
ies in the z ~ 1 sample with M{E) < -20, and the internal color dispersion of 
these galaxies simulated at z = 2.3. For each simulated galaxy, the /gu and 
//i6o images were constructed by resampling their V606 and /no images at a 
pixel scale corresponding to z = 2.3, and appropriately dimming their surface 
brightness. Note that the simulated C(/8i4,^i6o) may be underestimated be- 
cause this color at this redshift is slightly more sensitive to heterogeneity in a 
galaxy's stellar populations compared with ^(V^oe^Jiio) at z ~ 1, see §[" 



are occupied by local galaxies. However, it seems unavoid- 
able that the progenitors of the largest galaxies in the z ~ 1 
sample and at z ~ are not consistent with simple luminosity 
evolution. Pure luminosity evolution would fade the z ~ 2.3 
galaxies down the luminosity axis of Figure [3j where they 
have sizes that are too small compared with the observed dis- 
tribution of galaxy sizes at z ~ 1 . This result does not suffer 
from a selection bias in the sense that if the fully formed pre- 
cursors to the z ~ 1 galaxies were already in place at z = 2.3 , 
they would be observed in the deep NICMOS data (see § 13.21) . 
Therefore, the z ~ 2.3 population are not the fully-formed 
predecessors of the z ~ 1 galaxies. 

An alternative indication of growth for galaxies from z ~ 
2.3 to 1 is evident from their radius-stellar-mass distribu- 
tion, which is shown in Figure |9l Galaxy stellar masses have 
been estimated by fitting a suite of stellar-population syn- 
thesis models to the full spectrophotometry of the HDF-N 
galaxies (see Papovich et al. 2001; Dickinson et al. 2003). 
The stellar mass of a galaxy generally only increases with 
time (in the absence of interactions or tidal stripping that 
could eject stellar material),^ which implies that evolution in 
a galaxies' stellar population manifests itself only in the pos- 
itive direction of the ordinate in Figure |9l Galaxy-size evo- 
lution depends on the assembly history. Hierarchical models 
predict that the characteristic galaxy sizes scale roughly pro- 
portional to cosmic time (e.g., Bouwens & Silk 2002) from 
the accretion of surrounding material. This implies evolu- 
tion in the positive direction along the abscissa in Figure |9l 
by a roughly factor of 1.5 from z ~ 2.3 to 1. However, the 
final half-ma^^ size depends on the ratios of the sizes and 

^ As stars evolve off the main sequence they eject a portion of their stellar 
mass back into the galaxies' gas reservoir. This amounts to a decrease in the 
total stellar mass of ~ 30% over the lifetime of a passively evolving stellar 
population. 
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Fig. 9. — The galaxy size-stellar-mass distribution for the z ~ 1 and 2.3 
samples from the HDF-N. Stellar mass estimates correspond to the derived 
values from Papovich et al. (2001) and Dickinson et al. (2003). Error bars 
on the stellar masses are shown for one-third of the objects to illustrate char- 
acteristic uncertainties. The dashed horizontal line shows the stellar mass 
of a present-day -galaxy (see, Cole et al. 2001, with adjustments to the 
default cosmology). The triple-dot-dashed and dot-dashed lines show the 
limiting stellar-mass for a maximally old object observed at the (lOcr) flux 
limit of the data at z = 2.3 and 1, respectively. The data are sensitive to all 
sources with masses greater than these lines. Objects below these lines have 
lower mass-to-light ratios than the maximal model and are readily detected, 
but the completeness below these lines is less secure. 



masses of the progenitors. For example, using the relations 
presented in Cole et al. (2000) the radius of the merger prod- 
uct from two near-equal mass progenitors with size ratios, 
R2/R1 = (1/3, 1, 3), isRc^ (0.7, 1.3, 2.2) x Ru respectively. 
Note that in the case of R2/R1 = 1/3 that the final size is de- 
creased relative to the size of the largest progenitor. There- 
fore, in the case of major mergers evolution in either direction 
along the abscissa of the figure is possible. As illustrated in 
Figure |9| the distribution of z ~ 1 galaxies extends to both 
larger stellar masses and sizes than galaxies at z ~ 2.3. The 
larger-stellar-mass galaxies (M. ^ 5 x 10^^ A^©) at z ~ 1 are 
not present in the z ~ 2.3 sample. At the same time, galaxies 
with larger half-light radii (ri /2 ^ 5 kpc) appear at z ~ 1 in 
the HDF-N but not at higher redshifts. These galaxies in gen- 
eral have M < 10^^ A^0, similar to that found for the most 
massive galaxies at z ~ 2.3. The lack of larger galaxies with 
higher stellar masses at z ~ 1 is not unexpected simply be- 
cause such objects are very rare and unlik ely t o be present 
within the volume of the HDF-N (see also $ 14. 1ft . 

Previous studies of UV-bright z ~ 3 galaxies have 
concluded that their space density is roug hly that of 
present-day, ma ssive, early-type galaxies ( Steid el et al.ll996t 
iGiavalisco et al.l 1998). Through arguments of the measured 
z ^ 3 clustering and angular separation, this population may 
trace a one -to-one c orrespondence with massive dark matter 
halos (Adelberger et al. 1998: GiavaHsco & Dickinson 2001: 
iPorciani & Giavaliscoll2002h ~ Because the cores of such mas- 



DIVERSITY IN THE STELLAR-POPULATIONS OF GALAXIES 



13 



sive halos will presumably evolve to become the locations of 
massive present-day galaxies, it is natural to equate this high- 
redshift galaxy population with that of present-day "L*"- 
sized objects. One may expect galaxies in massive haloes to 
have an accelerated evolution, and thus they should assemble 
themselves at larger look-back times, which may partly be re- 
sponsible for some of the observed evolution. This scenario 
is plausible, with the additional requirement that the galax- 
ies necessarily undergo additional star formation and growth 
at z ^ 2 to match the z ~ 1 and present-day galaxy radius- 
luminosity and radius-stellar-mass distributions. 

4.2. The Distribution of Galaxy Colors from z ~ 1 - 3 

The scatter in the distribution of UV-optical total colors 
and internal color dispersion increases significantly in galax- 
ies from z ~ 2.3 to 1. FigurefTol shows color-color plots and 
the internal color dispersion as a function of total color for the 
z ^ 1 and 2.3 galaxy samples. The figure also shows the colors 
from models of an evolving stellar population at z = 1 and 2.3 
with three different monotonic star-formation histories (con- 
stant, instantaneous burst, and exponentially decaying) as a 
function of age from 10^-10^^ yr. 

The z ~ 2.3 galaxies span a relatively small range of to- 
tal color-color space even though they come from a IR- 
selected, flux-limited sample with complete photometric and 
spectroscopic redshift information (see § |2|l. The z ~ 2.3 
galaxies have colors that are broadly consistent with those 
of a young, constantly star-forming stellar population with 
moderate amounts of dust reddening (see also Sawicki & Yee 
[1998^, 'Papovich. Dickinson. & Ferguson 2001). Several ob- 
jects have highly reddened rest- frame UV colors (^450-/814 ~ 
1; indicated by open symbols in the figure). All of these 
sources have redshifts z ~ 2.9 at which point the flux in the 
5450-band is affected by Lyman a absorption from the IGM. 
This is not included in the model colors. Lyman a absorption 
does not affect the /g 14-/^^1 60 internal color dispersion. 

In contrast, the z ~ 1 galaxies show a greater scatter in 
their total colors. A significant fraction of the population 
has t/3oo-V606 and V606-'^iio colors that are well represented 
by stellar populations with declining star-formation histories, 
and ages greater than several gigayears. The z ~ 1 galaxies 
have colors spanned by the range of the models, including 
moderate dust extinction. Many of these objects have very red 
colors, ^/3oo- V6O6 > 1.5, V6O6-/110 > 2, and S/N([/3oo) < 10, 
which implies their emission at UV-optical wavelengths is 
dominated by older stellar populations with little contribution 
from recent star formation. Several galaxies have rest-frame 
UV colors that are apparently inconsistent with the models 
(i.e., galaxies with t/soo-^SH ^ 1-5 and Veoe-Jno < 1; in- 
dicated by open symbols in the plot). Two of these objects, 
HNM 748 and 1078, have spectroscopic redshifts 0.9, but 
photometric redshifts of 0.43 and 2.7, respectively. Thus the 
spectroscopic values (which were used to define the galaxy 
samples) are somewhat dubious. The other object, HNM 
1141, has a photometric redshift with maximum likelihood 
at Zph = 1-2, but with a poor goodness of fit. There is a sec- 
ond maximum near Zph ~ 3 that has a nearly-equal likelihood. 
Given their redshift ambiguities, we will exclude these objects 
in subsequent analysis. 

There is no clear correlation between total UV-optical color 
and internal color dispersion. While the ^(/8i4,^i6o) values 
of the z ~ 2.3 galaxies are generally low, the galaxies with 
highest internal color dispersion (^ > 0.05) have the reddest 
total UV-optical colors (/8i4-^i60 > 1 and5450-/8i4 > 0.5). 



This suggests that the red colors are due to either a mix of 
the stellar populations and/or variations in dust attenuation in 
the galaxies. At z ~ 1, many of the objects with high internal 
color dispersion have moderate rest-frame UV-optical col- 
ors, V606 -'^iio ^ 1- This is similar to the trend between the 
internal color dispersion and UV-optical color in local galax- 
ies (see P03, their figure 5), and suggests that high internal 
color dispersion results from a heterogeneous mix of young 
and old stellar populations (which produce the moderate to- 
tal UV-optical colors) that are spatially segregated within the 
galaxies. The increase in the internal color dispersion for the 
galaxies in the z ~ 1 HDF-N sample — as well as the in- 
creased scatter in total colors of the galaxy population — im- 
plies that the z ~ 1 galaxy population has more diversity in its 
stellar populations compared to z ~ 2.3. We interpret this as 
evidence for old stellar populations that are not prominent at 
z>2. 

Previous authors have studied the spatially resolved col- 
ors and color gradients of distant galaxies in order to inter- 
pret the star-formation histories of the stellar populations in 
these objects (e.g., Abraham etal. 1999; Corbinetal. 20Q2 
iMenanteau et al. 2001, 2004; Moth & Elston 2002). These 
studies have generally concluded that the spatial colors within 
galaxies reflect the evolutionary histories and internal dy- 
namics of the galaxies' stellar populations. Abraham e t al.l 
(1999) and Menanteau et al. (2001, 2004) argued that the in- 
ternal UV-optical colors are more sensitive to the relative 
ages of the galactic stellar populations than either dust or 
metallicity effects, and homogeneity in an object's internal 
color distribution is indicative of a broadly co-eval history 
for the galaxy's stellar populations. We reach a similar con- 
clusion in our analysis of the local galaxy sample (P03). 
A more diverse star-formation history for a galaxy's stellar 
populations produces heterogeneous internal colors. For ex- 
ample, by investigating the internal colors of high-redshift 
(z ^, 1) galaxies in HST images, both Abraham et al. ( 1999) 
and Menan teau et al.l (1200 1 ) find a large fraction (> 30-40%) 
of early-type galaxies in the field that show evidence for re- 
cent star formation activity. 

Clearly, galaxies with very strong internal color gradients 
and heterogeneities as studied by these authors would proba- 
bly also have large internal color dispersion. A small amount 
of ongoing star formation in an already massive galaxy can 
dominate its UV flux emission and give rise to strong varia- 
tions in a galaxy's internal colors, even though it may con- 
stitute a negligible total to the galactic stellar mass (e.g., the 
"frosting" scenario of young stars in elliptical and lenticular 
galaxies; see Trager et al. 2000). Our measure of the in- 
ternal color dispersion has an advantage that it can quantify 
the amount of ongoing star forming relative to the previously 
formed populations (see discussion below). 

4.3. Simulations of Galaxy Internal Colors 

A fraction of the z ~ 1 galaxies in HDF-N have morpholo- 
gies of early-to-mid-type spiral galaxies, and these objects 
exhibit the highest internal color dispersion (see figures |6|l. 
This agrees with trends between morphology and the internal 
color dispersion observed in local galaxies (P03). 

In most hierarchical merging models, the growth of ro- 
tationally supported galaxy disks is expected as gas cools 
and accretes towar d the center of th e underlying dark- 
matter halos (e.g.. Mo. Mao . & White 1998), and this is 
broadly supported by the evolution in galaxy sizes (see 
§ 14.111 . This process is expected to form star-forming disks 
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Fig. 10. — Distribution of galaxy internal color dispersion and total colors of galaxies at z ~ 1 (top panels) and z ~ 2.3 (bottom panels) in the HDF-N. The 
left panels show the derived internal color dispersion between the UV-optical bands as a function of the total UV-optical color. The blue, dashed horizontal lines 
indicate a fiducial value ^ = 0.05. The right panels show the (total) color-color plots for each galaxy sample. The open stars in the z ~ 1 plot correspond to objects 
with insecure redshifts (see text). The open boxes in the z ~ 2.3 plot correspond to objects with z ~ 2.9, where the 5450-band flux is attenuated by intergalactic 
Lyman a, which increases the observed ^450-/8 14 color. In these plots, the three curves illustrate the colors in the observed frame for an evolving stellar 
population with three monotonic star-formation histories in the observed frame for z = 1 (f/soo ^606-^1 10) and z = 2.3 (5450/814/^^160) respectively. The blue curve 
with star symbols corresponds to a stellar population with a constant star formation rate. The green curve with triangle symbols corresponds to a stellar population 
with an exponentially decaying star-formation rate and e-folding time constant, r = 1 Gyr. The red curve with circles corresponds to a stellar population formed 
in an instantaneous burst. The symbols along each curve indicate the colors of the stellar population at an age of log(f/yr) = 7,8,8.5,8.7,9,9.5,9.7, 10. The 
arrows illustrate the reddening vector for a color excess, E(B-V) = 0.2, and assuming a starburst-like extinction law (Calzetti et al. 2000). 



around spheroids produced in merger remnants, and cre- 
ate spiral galaxies, which broadly reproduces the observed 
distribution of colors, gradients, and masses in z ^ 2- 
3 ga l axies dC ontardo. Steinmetz. & Frit ze-von AlvenslebenI 
ll998t 'Steinmetz & Navarro 2002). 

To investigate how the formation of new disks around older 
spheroids affects the internal color dispersion, we considered 
a simple model to represent the galaxies at z = 1 and 2.3. We 
simulated galaxies with two components, an old spheroidal 
component and a young star-forming disk component, and 
assigned to these colors using the Bruzual & Chariot (2003) 
stellar population synthesis models and assuming a Salpeter 



IMF and solar metallicity (see Table [TJ. The spheroidal com- 
ponent was given the colors of a stellar population formed 
in an instantaneous burst viewed at an age of 2.8 Gyr and 
a r^/^-law surface-brightness profile. The selected age cor- 
responds to the maximum possible at z = 2.3, and approxi- 
mately to the elapsed cosmic time between z = 23 and 1, and 
the half-light radius is consistent with the observed proper- 
ties of the z ~ 2.3 galaxy sample and with the derived sizes 
of local galaxy bulges (Carollo 1999). The disk component 
has the colors of relatively young stellar populations (see Ta- 
ble [Q, and an exponential surface-brightness profile. The 
disk brightness is incrementally increased relative to that of 
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Fig. n . — Distribution of galaxies internal color dispersion as a function of total UV-optical color. Data and symbols definitions are the same as in FigureFTol 
The curves in the correspond to the simple models discussed in the text. The black curves show the expected colors of galaxies with an old central bulge and 
star-forming exponential disk with different star-formation histories (model 1, dashed lines; model 2, solid lines; see text). The red curves illustrate the expected 
colors for objects with additional star-forming, H II regions within older, passively evolving stellar populations in a spheroid and disk (model 3). 



the spheroid to span the full range of bulge-to-disk ratios. 

The total colors and internal color dispersion of these mod- 
els are plotted in Figure [H] with the data. Interestingly, the 
same model predicts larger ^(/8i4,//i6o) values at z = 2.3 rel- 
ative to ^(V606,^iio) at z = 1 by ^ 40%. The ku-Hm color 
at z = 2.3 is bluer in the rest frame and covers a longer rest- 
frame wavelength baseline than the Veoe-Jno color at z = 1. 
Therefore, the I^u -//i6o is slightly more sensitive to the UV- 
optical color dispersion, and an object with the same intrinsic 
color dispersion would have a greater measured ^(/8i4,//i6o) 
3.tz = 2.3 than ^(V^oe^Jno) at z = 1. This actually strengthens 
the conclusion that the range of the UV-optical internal color 
dispersion increases from z ~ 2.3 to 1. 

The simple models of a star-forming disk surrounding an 
older central bulge span the range of total color and internal 
color dispersion observed in most of the galaxies at z ~ 1 that 
have moderate UV-optical colors. At low and high disk-to- 
bulge ratios, the model colors are nearly uniform (i.e., domi- 
nated by the disk and bulge respectively). The internal color 
dispersion is largest for mixed models. The stellar populations 
in these galaxies are consistent with this scenario. 

However, there are many galaxies in the z ~ 1 sample with 
very red colors (V606 --^iio ^ 1-5) and high internal color dis- 
persion, and these are not accounted for by the simple mod- 
els. The light from these objects must be dominated by rela- 
tively evolved stellar populations, possibly with small traces 
of ongoing star-formati on to in duce large internal UV-optical 
color dispersion (Trag er. et al.l l200Q). To model this possi- 
bility, we considered the spheroid and disk models with an 
additional source of very young star-forming knots ("H II re- 
gions") typically located in galactic spiral arms and starbursts 
(described in Table [I}. These models, plotted in Figure fTTl 
reproduce the high internal color dispersion with redder to- 
tal colors than the pure disk/bulge models. The high inter- 
nal color dispersion is produced from spatially segregated and 
heterogeneous stellar populations. We therefore conclude that 
the increase in the internal color dispersion between z ~ 2.3 
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Fig. 12. — The measured internal UV-optical color dispersion from simu- 
lated galaxies as a function of the fraction of stellar mass in the young stellar 
populations of the disk, i.e., disk / disk + bulge)- The three curves cor- 
respond to galaxies with the spheroidal components and disk components 
listed in Table [Tl (model 1, dashed line; model 2, thick line; model 3, thin 
line) viewed at z = 1 . 



to 1 corresponds to greater diversity in the amounts of young 
and old stellar populations in these galaxies. 

The best example of this behavior in the z ~ 1 HDF-N sam- 
ple is the giant, "grand-design" spiral galaxy, HNM 1488 (see 
figures [T] and 01. Bulge and possible bar features dominate 
the rest-frame 5-band (/no image). In the rest-frame UV 
(V606), these features fade relative to the disk, spiral arms, and 
prominent star-forming knots embedded in the disk. The red- 
der bulge dominates the rest- frame optical light, giving a red 
total color, V606 --^iio =1.8. The internal color dispersion is 
very high, (,(V^o6^Jno) = 0.13, and arises primarily from the 



16 



PAPOVICH ET AL. 



TABLE 1 

Galaxy Component Models 









ri/2 




M/Lb 


^606 --^110 


/814-^160 


Component 


t 


t/r 


[kpc] 


Profile 


[Mq/Lb,q] 


z=l 


z = 2.3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


disk model 1 


0.3 Gyr 


0.6 


7.5 


Exponential 


0.11 


0.41 


0.36 


model 2 


0.1 Gyr 


1.0 


7.5 


Exponential 


0.065 


0.15 


0.08 


model 3 


0.1 Gyr 


4.0 


7.5 


Exponential 


0.26 


1.45 


1.67 


spheroid 


2.8 Gyr 


oo 


2.3 


de Vaucouleurs 


2.7 


3.06 


4.83 


H II region 


50 Myr 


^0 




Point source 


0.032 


-0.16 


-0.27 



Note. — (1) Galaxy component; (2) Age of stellar population; (3) Ratio of age to the 
star-formation rate e-folding timescale; (4) Half-light radius; (5) Profile; (6) Stellar mass 
to blue light ratio in solar units; observed frame colors at (7) z = 1, and (8) z = 2.3. 
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Fig. 13. — Distribution of concentration and asymmetry parameters of the z ~ 1 (Left panel) and z ~ 2.3 (Right panel) galaxy samples. Filled stars denote 
objects with ^(Veoe^Jno) > 0.05 and ^(/8i4,//i6o) > 0-05 in the z ~ 1 and z ~ 2.3 panels, respectively. Objects with internal color dispersion values below this 
value are indicated by filled circles. The dashed lines denote rough fiducial regions in the C-A diagram that correspond to early (ellipticals, early-type spirals), 
middle (mid-type spirals), and late (irregulars and late-type spirals) Hubble-type galaxies (see Conselice et al. 2000). 



color differences between the bulge, disk, and H II regions 
(Figure 0. Based on the models, these three features are re- 
quired to produce both the red total color and high internal 
color dispersion. 

The different galaxy components described in the simula- 
tions above have very different stellar-mass-to-blue-light ra- 
tios, as indicated in Table [l] Figure [T2| shows the internal 
color dispersion of the models as a function of the relative 
amount of stellar mass in the different model components. 
For the disk models listed in Table [I] the stellar mass of the 
bulge exceeds that in the disk by up to a factor of ~ 10-40 
when the luminosities of the bulge and disk are equal, and the 
internal color dispersion reaches a maximum near this point. 
Based on Figure [121 high values of the internal color disper- 
sion > 0.05) occur when < 10% of the stellar mass resides 
in recently formed stars (see also P03). 

4.4. Internal Color Dispersion as a Function of Morphology 

The increase in range of the UV-optical internal color 
dispersion coincides with the emergence of large, early- 

and mid- type Hubble-seq uence galaxie s. For example, 

[Abraham et al. ( 1999) and Conselice et al.| (|2003|) have shown 
that galaxy morphological type is generally distinguished 



by their co ncentration and asymmetry indexes (see also 
IConselicJ2003J . If the formation of Hubble-sequence galax- 
ies corresponds to an increase in the scatter of the UV-optical 
color dispersion, then one expects there to be a relation be- 
tween the internal color dispersion and these morphological 
indicators. FigureElshows the concentration and asymmetry 
values along with int ernal co lor dispersion for the two HDF- 
N samples (Conseli ce et al.l l2003). At z ~ 2.3 galaxies with 
high UV-optical internal color dispersion occupy the upper 
range of asymmetry values, which is indicative of irregular 
star-forming galaxies possibly due to interactions and merg- 
ers (Conselice. Blackburne. & Papovich 2004a). This sug- 
gests that the internal color dispersion results from variations 
in the stellar populations and dust opacity, e.g., HNM 813 + 
814 (see figures |2| and 0. Because galaxies with relatively 
high internal color dispersion correspond to highly asymmet- 
ric galaxies, this supports the notion that they are strongly 
interacting and/or merging. At z ~ 1 galaxies with high inter- 
nal color dispersion span a range of morphological types as 
traced by concentration and asymmetry. There is still a group 
of galaxies with high internal color dispersion and high asym- 
metry. These objects are similar to the higher redshift star- 
forming galaxies, and is consistent with a recent analysis of 
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HST/ ACS images, which suggests that star-forming galaxies 
at z ~ 1.5 and 4 have similar far-UV (1500 A) morphologies 
(Lotz et al. 2005). There are also a large number of galaxies 
with low asymmetry, and higher concentration, which suggest 
that the high internal color dispersion in these objects occurs 
in otherwise normal mid/early-type galaxies. It seems that 
the emergence of Hubble-sequence galaxies coincides with 
larger internal color dispersion, and thus the assembly of di- 
verse, spatially segregated stellar populations. 

Two of the galaxies, HNM 1214 and 1453, are classi- 
fied as early-types by high concentration and low asymme- 
try, and also have moderately high internal color dispersion 
:^ 0.08). Upon visual inspection in both cases very faint 
blue companions are visible at distances of < V\ which were 
not split by SExtractor and increase the internal color disper- 
sion. The companion near HNM 1214 is visible roughly one 
arcsecond to the left of the galaxy in Figure [I] Companions 
only affect the internal color dispersion in very few of ob- 
jects (these two and HNM 145 3, wh ich was eliminated from 
the sample for this reason, see $ 13.3ft . The large internal color 
dispersion in the majority of objects results from spatially seg- 
regated heterogeneous stellar populations. 

4.5. Formation of Stellar Diversity in Galaxies from z ~ 1 - 3 

Many of the z ~ 1 galaxies have red total colors that are con- 
sistent with the passively evolving stellar-population models 
with ages older than 1-5 Gyr. Because these ages are greater 
than the elapsed cosmic time from z ~ 2.3 to 1, some of the 
stellar populations currently in z ~ 1 galaxies must have been 
present in progenitors at the earlier epoch. However, due to 
the large growth in galaxy stellar mass and size over these red- 
shifts, the z ~ 1 galaxies were generally not fully formed by 
z>2. 

In addition to evolution in size and mass, the increase in 
the range of galaxy internal color dispersion from z ^ 23 
to 1 is evidence that the diversity in the stellar populations 
in galaxies is increasing. Papovich et al. (2001) found that 
the relatively short ages of galaxies at 2 < z < 3.5 (less than 
a few hundred megayears) coupled with the lack of older 
faded "post-burst" remnants implied that these galaxies have 
recurring star-formation episodes. The time scales of these 
episodes correspond to < 1 Gyr such that the newly-formed 
stellar populations dominate the rest-frame UV-optical wave- 
length regions of the galaxy SEDs. One cannot exclude the 
possibility that older (and possibly more heterogeneous) stel- 
lar populations are hidden beneath the "glare" of the young 
stars. This is broadly consistent with small range of low inter- 
nal color dispersion observed in the z ~ 2.3 sample, but with 
the additional constraint that the nascent stellar populations 
in these galaxies dominate both the total galaxy colors and 
the internal galaxy colors, at least over scales of :^ 1.5 kpc 
(the NICMOS resolution at these redshifts). 

If either recurrent bursts such as those described above 
or quasi-constant star formation is the norm for galaxies at 
z ~ 2.3, then eventually these galaxies will form sufficient 
stellar mass, which will dominate the emitted optical lumi- 
nosity. Strong UV-optical internal color dispersion generally 
requires that > 80% of t he to tal stellar mass resides in older 
stellar populations (see $ l4.3li . At z = 2.3, the age of the Uni- 
verse is less than 3 Gyr, and thus there has been little opportu- 
nity for enough star-formation events to have instilled strong 
heterogeneity in the galaxies' stellar populations. At z ~ 1 
the age of the Universe roughly doubles to ~ 6 Gyr, which 
permits a greater number of star-formation episodes with this 



timescale and allows for greater internal color dispersion. The 
increase in the internal color dispersion for z ~ 2.3 to 1 seems 
to require that enough time has elapsed for galaxies to have 
formed a sufficient amount of their stellar mass relative to 
their current star-formation rate such that the older stars dom- 
inate the optical luminosity. 

The discussion above ignores the effects of dust extinction 
in star-forming regions, which can contribute to the over- 
all "patchiness" in galaxy inter nal colors and produc e strong 
morphological ^-corrections (Windhorst et al. 2002 ). How- 
ever, studies of the variation in internal colors ( Abrahamet alj 
1999; Menanteau et al. 2004) and the internal color dispersion 
of local galaxies (P03) broadly conclude that these diagnos- 
tics are more sensitive to variations in the ages of the galaxy 
stellar populations relative to variable dust opacity. There- 
fore, we attribute the strong variation in galaxy internal color 
dispersion primarily to differences in the ages of the galaxy 
stellar populations. 

Bec ause the half-light radius is growing from z ~ 2.3 to 1 
(see § 14.111 . we know that galaxies do not form in situ and 
evolve passively over this redshift range. Furthermore, as the 
size increases, galaxies at later epochs will have proportion- 
ally more young stars at larger radii. So, size evolution will 
naturally lead to color gradients. This is supported by the fact 
that the z ~ 1 galaxies with the largest internal color disper- 
sion appear to have red centrally concentrated cores with bluer 
regions at larger radii (see Figure |6ll. 

Although the large internal color dispersion values at z ~ 1 
imply a large stellar diversity at this epoch, the same is not 
t he ca se at z ~ 2.3. Based on the simple models presented in 
$ 14.31 although there is sufficient time at z = 2.3 for galaxies 
to have developed a heterogeneous stellar content with sig- 
nificant internal color dispersion, none is observed (see, for 
example. Figure [HJ. Therefore, we conclude that galaxies at 
z ^ 2 are not forming a diverse stellar content as rapidly as is 
allowed purely from timescale arguments. This implies that 
galaxies at these epochs are not assembling their stellar pop- 
ulations in the quasi-continuous, inside-out models described 
in §10 

One alternative to account for the fact that the z ~ 2.3 
galaxies do not exhibit strong internal color dispersion is 
that intense, discrete starbursts such as those associated 
with actively merging galaxies dominate the star-formation 
histories at this epoch. Intense starbursts naturally dom- 
inate the UV-optical light at the time the galaxy is ob- 
served, which accou nts for the obse rved prope rties of the 
galaxie s at z ^ 2 - 4 ("Lowenthal et al. 1997; Sawicki & Yee 
J 998: iSt eidel et al.| [l ^99: Papovich. Dickinson. & Ferguson 
l200lUPapovich et al.ll2004 . Furthermore, if these starbursts 
arise from strong interactions or mergers, then during the 
merger one expects some internal color dispersion to re- 
sult from any previously formed, heterogeneous stellar pop- 
ulations with some contribution from varying dust opacity 
(see, e.g., the large heterogeneity in both th e stellar popula- 
tions and dust extinction in NGC 403 8/4039 i lwhitmore et"an 
Il999h . The merger process will tend to erase any diver- 
sity in the constituent stellar-populations by fully mixing 
them on short timescales (e.g., Mihos & Hernquist 1996; 
Steinmetz & Navarroj l2002h . and thus the merger product 
should have largely homogenized stellar populations. This is 
consistent with the small derived values for the internal color 
dispersion at z ~ 2.3. The period while galaxies are strongly 
interacting and/or actively merging may produce patchy star 
formation with larger internal color dispersion. Indeed, the 
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objects at z ~ 2.3 with high internal color dispersion ~ 
0.05 -0. 1) have large asymmetries (FigurefT3l. supporting the 
hypothesis that they are actively merging. We therefore inter- 
pret the observations as evidence for major mergers associated 
with strong starbursts. Furthermore, these events probably oc- 
cur frequently, with the timescale between bursts < 1 Gyr, in 
order to keep the stellar content in these galaxies homoge- 
nized and to keep the total UV-optical colors dominated by 
the young starburst. 

By z ~ 1 the internal color dispersion in galaxies has in- 
creased significantly, which we interpret as an increase in the 
diversity of galaxies' stellar content. This occurs partly be- 
cause galaxies have formed sufficient mass in older stellar 
populations that recent star forming episodes contribute less 
to the total optical emission (while still dominating the UV 
light; see the discussion above). However, the heterogeneous 
stellar populations are spatially segregated: the z ~ 1 galaxies 
with the highest internal color dispersion typically have red 
central cores surrounded by bluer regions (see Figure|6|). This 
signifies that major mergers occur less frequently, because any 
merger activity would tend to erase this diversity in the con- 
figuration of the stellar populations. The growth of large disks 
at this epoch naturally increases the internal color dispersion, 
and is consistent with the interpre tation of the ubiquitqusnes s 
of thick disks in local galaxies ( Dalcanton & Bernstein'2002V 
Thus, the emergence of galaxies with large internal color dis- 
persion signifies great diversity in their stellar content, and we 
interpret this as evidence for a decline in the major-merger 
rate in galaxies at z 5, 1-5 (see als o van den Bergh et al. 2002; 
IConselice. Blac kburne. & Papovich 2004a). 

The decline in the importance of bursts with decreasing red- 
shift may imply that more quiescent modes dominate the total 
star-formation rate. This is consistent with the suggestion by 
van den Bergh (2002) and reminiscent of the "collisional star- 
burst" models of Somerville et al. (2001), where bursts domi- 
nate at high redshifts, and quiescent star formation dominates 
at lower redshifts (z ^ 0.8). The larger range of internal color 
dispersion for the massive galaxies at z ~ 1 implies that quies- 
cent star-formation plays a more dominant role at later times 
compared to major mergers. The transition between where 
these modes dominate occurs between 1 .4 - 1 .9. 

5. CONCLUSIONS 

We have studied the evolution in the morphologies and col- 
ors of two samples of high-redshift galaxies at z ~ 2.3 and 
z ^ 1 in order to investigate how these galaxies assembled 
their stellar populations. From z ~ 2.3 to 1, galaxies grow 
in size, stellar mass, total color, and internal color dispersion, 
which we interpret as evidence for an increase in the diversity 
of stellar populations as galaxies assemble. 

The majority of luminous galaxies (M(B) < -20.5 in the 
rest frame) in the z ~ 1 sample have rest-frame optical mor- 
phologies that are classifiable as early-to-mid Hubble types. 
These galaxies have regular, symmetric morphologies, and 
many show strong transformations between their morpholo- 
gies when viewed at rest-frame UV and optical wavelengths. 
In contrast, the rest-frame optical morphologies of the lu- 
minous (M(B) < 20.5 in the rest-frame) galaxies at z ^ 23 
are irregular and/or compact, and there is little difference in 
the galaxy morphologies from rest-frame far-UV to optical 
wavelengths. None of the z ~ 2.3 galaxies appear to be nor- 
mal Hubble-sequence galaxies. Because the z ~ 2.3 galaxies 
show little transformation in morphology from rest-frame UV 
to optical wavelengths, the UV components are not generally 



the small, star-forming pieces within larger host galaxies. 

The mean galaxy size increases from z ~ 2.3 to 1 by roughly 
40%, in broad agreement with expectations from hierarchi- 
cal models. Furthermore, the number density of large galax- 
ies, ri/2 > 3 kpc and M(B) < -20, increases by a factor of 
^ 7. We have tested that the size evolution is robust against 
surface-brightness dimming effects. The half-light radii of 
the z ~ 1 galaxies when simulated at z = 2.7 decrease by 
< 20% compared to their measured values at z ~ 1 , but this 
decrease is small compared to the observed evolution. The 
radius-luminosity distribution of galaxies at z ~ 1 is consis- 
tent with the local distribution for pure luminosity evolution. 
However, the size-luminosity, and size-mass distributions of 
z ~ 2.3 galaxies indicate that they are not the fully formed 
progenitors of large, present-day galaxies. Galaxies at z ^ 2 
must continue to build, both in terms of size and stellar mass. 

We have analyzed the rest-frame UV-optical total colors, 
and the internal color dispersion using the novel statistic de- 
veloped in P03. The internal color dispersion quantifies dif- 
ferences in galaxy morphology as a function of wavelength, 
and between rest-frame UV and optical wavelengths it con- 
strains the amount of current star-formation relative to the 
older stellar populations. Both the mean and scatter of total 
color and the internal color dispersion increase from z ~ 2.3 
to 1 . At z ~ 1 , galaxies with high internal color span a range 
of morphological types. Many are spiral galaxies, with a few 
starbursting and interacting systems. There is no clear cor- 
relation between total color and the internal color dispersion, 
although galaxies with the highest internal color dispersion 
have moderate total UV-optical colors, which implies a mix 
of heterogeneous stellar populations. The z ~ 1 galaxies with 
the highest internal color dispersion appear to be early-to- 
mid-type spiral galaxies, and they exhibit strong transforma- 
tions between their UV and optical rest-frame morphologies. 
At z ~ 2.3 few galaxies have high internal color dispersion, 
and those that do also have interacting and disturbed mor- 
phologies, which implies that the UV-optical light from these 
galaxies is dominated by young, largely homogeneous stellar 
populations. We interpret the evolution in galaxy color and 
internal color dispersion as evidence that the diversity in the 
stellar populations of galaxies is increasing from z ~ 2.3 to 1 
as the older stellar components build up over time. 

To interpret the evolution in the properties of galaxies, we 
have modeled the stellar populations of galaxies in spheroidal 
and disk components, and star-forming H II regions. These 
simple models broadly support the conclusion that the range 
of internal color dispersion corresponds to increased diver- 
sity in the galaxies' stellar populations. Galaxies with signif- 
icant bulge components with star-forming disks and/or com- 
pact, H II regions produce high internal color dispersion, and 
broadly reproduce the colors, sizes, and internal color disper- 
sion of the high-redshift galaxies. We find that large values 
of the internal color dispersion require galaxy to have formed 
spatially segregated, diverse stellar populations where old 
stars dominate the optical light and stellar mass (i.e., > 80% 
of the stellar mass resides in stars > 1 Gyr). For smaller frac- 
tions of old stars, young stellar populations dominate the UV- 
optical internal colors. 

The scatter in the internal color dispersion of the z ~ 2.3 
galaxies is smaller than what is allowed under these models 
and basic timescale arguments. We interpret this as evidence 
that brief, discrete and recurrent starburst episodes dominate 
the star-formation history of galaxies at this epoch. If these 
arise from strong interactions or mergers of gas-rich con- 
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stituents, then they will erase any heterogeneity in the stel- 
lar content that otherwise develops, and this is supported by 
the low internal color dispersion observed in galaxies at this 
epoch. The greater range of internal color dispersion that 
forms at z ~ 1 suggests that major mergers are less frequent 
and more quiescent star-formation mechanisms are the norm. 

The greater diversity in the stellar populations of these 
high-redshift galaxies coincides with the emergence with 
large Hubble sequence galaxies. Large values of the UV- 
optical internal color dispersion require a diverse and spatially 
unmixed stellar population, and this occurs when galaxies 
have formed most of their stellar mass. In order to maintain 
the spatial heterogeneity between the young and old stellar 
populations also requires that major mergers are less common 
at this epoch. These conditions all occur at z ^ 1 .4, which nat- 
urally allows galaxies with Hubble-sequence morphologies to 
develop. 
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